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ABSTRACT 

We use two-dimensional hydrodynamic simulations of viscous disks to examine whether dynamically- 
interacting multiple giant planets can explain the large gaps (spanning over one order of magnitude 
in radius) inferred for the transitional and pre-transitional disks around T Tauri stars. In the absence 
of inner disk dust depletion, we find that it requires three to four giant planets to open up large 
enough gaps to be consistent with inferences from spectral energy distributions, because the gap 
width is limited by the tendency of the planets to be driven together into 2:1 resonances. With 
very strong tidal torques and/or rapid planetary accretion, fewer planets can also generate a large 
cavity interior to the locally formed gap(s) by preventing outer disk material from moving in. In 
these cases, however, the reduction of surface density produces a corresponding reduction in the inner 
disk accretion rate onto the star; this makes it difficult to explain the observed accretion rates of 
the pre/transitional disks. We find that even with four planets in disks, additional substantial dust 
depletion is required to explain observed disk gaps/holes. Substantial dust settling and growth, with 
consequent significant reductions in optical depths, is inferred for typical T Tauri disks in any case, 
and an earlier history of dust growth is consistent with the hypothesis that pre/transitional disks are 
explained by the presence of giant planets. We conclude that the depths and widths of gaps, and 
disk accretion rates in pre/transitional disks cannot be reproduced by a planet-induced gap opening 
scenario alone. Significant dust depletion is also required within the gaps/holes. Order of magnitude 
estimates suggest the mass of small dust particles (< l/im) relative to the gas must be depleted to 10~^ 
- 10~^ of the interstellar medium value, implying a very efficient mechanism of small dust removal or 
dust growth. 

Subject headings: accretion disks, stars: formation, stars: pre-main sequence 
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1. INTRODUCTION 

The transitional and pre-transitional disks around 
young stars exhibit strong dust emission at wavelengths 
> 10/im, while showing significantly reduced fluxes rela- 
tive to typical T Tauri disks at shorter wavelengths (e.g., 
Calvet et al. 2002, 2005; D'Alessio et al. 2005; Espaillat 
et al. 2007, 2008). In the pre-transitional disks, there 
is evidence for emission from warm, optically-thick dust 
near the star (Espaillat et al. 2007, 2008, 2010), while in 
the transitional disks the emission at < 10/im appears to 
be due entirely to optically-thin dust (Calvet et al. 2002, 
2005; Espafllat et al. 2010). The depletion of near- to 
mid-infrared emission is generally interpreted as being 
due to evacuation of the disk interior to scales 5 to 
~ 50 AU (Marsh & Mahoney 1992; Calvet et al. 2002, 
2005; Rice et al. 2003; Schneider et al. 2003; EspaiUat et 
al. 2007, 2008, 2010; Hughes et al. 2009), an interpreta- 
tion confirmed in some cases via direct sub-mm imaging 
(e.g., Pietu et al. 2006; Brown et al. 2007, 2009; Hughes 
et al. 2009; Andrews et al. 2009). 

One proposed mechanism for clearing inner disks while 
leaving the outer disk relatively undisturbed is the forma- 
tion of giant planets, which can open gaps in disks (e.g.. 
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Lin & Papaloizou 1986; Marsh & Mahoney 1992; Nelson 
et al. 2000; Calvet et al. 2002; Rice et al. 2003), and 
which are expected to form initially at relatively small 
disk radii due to the shorter evolutionary timescales 
compared with the outermost disk. However, there are 
three observational challenges, as discussed further in 
§2, which theoretical explanations of these systems must 
confront. 

First, the transitional and pre-transitional disk systems 
exhibit average gas accretion rates close to T Tauri disk 
accretion rates (~ lO~^M0yr^^; Hartmann et al. 1998) 
onto their central stars (e.g. Calvet et al. 2002, 2005; 
Espafllat et al. 2007, 2008; Najita et al. 2007). Main- 
taining this accretion requires either a significant mass 
reservoir interior to the disk-clearing planets, or some 
way of allowing mass from the outer disk to move past 
the gap-clearing planets. 

Second, the cleared regions in these disks are large (e.g. 
Espaillat et al. 2010). In the case of the transitional 
disks, the optically-thin region must extend from radii 
as large as tens of AU all the way in to the central star. 
Even the pre-transitional disks, which have evidence for 
optically-thick dust emission in the innermost regions, 
must have large disk gaps. Furthermore, the spectral 
energy distribution (SED) modeling suggests that the 
detected optically thick region may only extend to radii 
~ 1 AU, with an extremely dust free region beyond, until 
reaching the outer optically-thick disk (Espaillat et al. 
2010; §2). 

Third, the requirement that the gap/hole be optically 
thin implies that the mass of dust in sizes of order a 
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micron or less must be extremely small. Thus, either the 
planet-induced gap is very deep and is effectively cleared 
of gas and dust, or the dust abundance is reduced by 
many orders of magnitude from abundances in the diffuse 
interstellar medium. 

These three conditions must be fulfilled simultane- 
ously. Photoevaporation has been proposed as one 
gap/hole clearing mechanism (Alexander & Armitage 
2007, 2009) , but requires the mass loss rate by the photo- 
evaporation to be comparable to the disk accretion rate. 
Initial estimates of photoevaporative mass loss suggested 
values of - 4xlO-i"M0 yr'^ (Clarke et al. 2001), which 
would not be rapid enough to counteract accretion. More 
recent estimates suggest higher values, perhaps as much 
as ~1O~^M0 yr~^ due to the inclusion of X-rays (Gorti & 
Hollenbach 2009; Owen et al. 2010); however, the highest 
values are problematic, because then it becomes difficult 
to understand why so many T Tauri disks last for sev- 
eral Myr with accretion rates smaller than 1O~*M0 yr~^ 
(Hartmann et al. 1998). 

Najita, Strom, & Muzerolle (2007) suggested that gi- 
ant planets might explain the transitional disks by cre- 
ating a gap while still maintaining accretion onto the 
central star. Najita et al. argued that there is some ev- 
idence of reduction in the accretion rates in transitional 
disks relative to the so-called "primordial" disks, per- 
haps by an order of magnitude, and that this was con- 
sistent with simulations of a single giant planet in the 
studies by Lubow & D'Angelo (2006) and Varniere et al. 
(2006). However, as discussed in §2, the accretion rates 
of the pre/transitional disks are not very low in absolute 
magnitude. Moreover, while Varniere et al. (2006) were 
able to produce a large cleared inner region with a single 
planet, the study of Crida & Morbidelli (2007) calls this 
result into question (see also §3). 

In this paper we consider whether disk gap formation 
by a multiple planet system can satisfy the following re- 
quirements: 1) creation of a gap extending over a large 
range in radius; 2) maintenance of the inner disk accre- 
tion rate by flow through the planet-created gap; and 
3) sufficient reduction in the disk surface density within 
the gap such that extreme depletion of disk dust is not 
essential. We examine the problem in the context of vis- 
cous disks, with improvements over some previous simu- 
lations, including a more realistic temperature distribu- 
tion and an improved inner boundary condition. We find 
that while multiple giant planets can indeed open large 
gaps, it is difficult to explain the inferred properties of 
the pre/transitional disks without invoking substantial 
depletion of small dust. We offer a speculative scenario 
in which combined dust depletion and planet formation 
explains the observations. 

2. INFERRED PROPERTIES OF PRE/TRANSITIONAL 
DISKS 

We first present a summary of properties of some of 
the known pre/transitional systems in Table 1 as derived 
from analysis of optical, infrared, and millimeter SEDs, 
including the pre-transitional disk systems LkCa 15, UX 
Tau A, Rox 44, and of the transitional disk GM Aur. 
The parameters are mostly taken from Espaillat et al. 
(2010), where the methods used to derive these estimates 
are discussed. 

As indicated in Table 1, the estimated maximum 



extent of the inner optically-thick region in the pre- 
transitional disks is ^ 0.2 — 0.4 AU, while the inner edges 
of the optically-thick outer disks lie at ~ 40 — 70 AU. The 
location of the inner edge of the outer disk inferred from 
SED modeling is consistent with that derived using mil- 
limeter interferometry (Pietu et al. 2006; Andrews et al. 
2009; S. Andrews, 2010, private communication). The 
maximum extent of the inner disk is less certain and is 
based on SED modeling of the near-infrared emission. In 
this wavelength region, the flux in pre-transitional disks 
is dominated by the emission of the inner wall. To limit 
the size of the optically thick disk behind the inner wall, 
Espaillat et al. (2010) varied the maximum extent of this 
inner disk and fitted it to the SED of the object. The 
value listed in Table 1 represents the largest Rinnerdisk 
that does not produce too much excess emission, assum- 
ing a typical disk flaring structure. A caveat is that if 
the inner wall is "puffed up" relative to the disk behind 
it (e.g. Natta et al. 2001), or if the dust behind the 
inner wall is more settled than the wall, the inner disk 
would be in shadow and not emit strongly. However, we 
note that the near-infrared interferometric observations 
by Pott et al. (2010) are consistent with the dust emis- 
sion being confined to a small radius as inferred from the 
SED modeling. The amount of emitting optically thin 
sub-micron-sized dust is small, ^ IO^^^Mq. 

Transitional disks do not have the inner optically thick 
disk observed in the pre-transitional disks. The tran- 
sitional disk of GM Aur has a 20 AU hole, based on 
modeling both the SED and submillimeter interferomet- 
ric visibilities (Hughes et al. 2009). About lO^^^ Mq of 
sub-micron-sized optically thin dust is necessary within 
the inner AU of the hole to reproduce the observed sil- 
icate emission in the near-infrared. Calvet et al. (2005) 
inferred that the inner emitting, optically thin dust in 
GM Aur is confined to within 5 AU by modeling its 
SED. While SED modeling is nonunique, the interfer- 
ometric observations of Pott et al. (2010) indicated that 
this dust is confined within even smaller radii, < 0.15 AU 
(see also Akeson et al. 2005). Thus, as stated by Pott 
et al. (2010), the pre/transitional disks are often charac- 
terized by distinct dust zones "which are not smoothly 
connected by a continuous distribution of optically thick 
material", which may be interpreted as the result of a 
highly cleared gap with an inner region of optically thin 
or thick dust. 

The final important point to be made is that, as shown 

in Table 1, the gas accretion rates onto the central stars 
are substantial, ^ 0.3 — 1 x IO^^Mq yr~^, not much lower 
than rates typical of T Tauri stars (e.g., Hartmann et al. 
1998). These high gas accretion rates agree with signif- 
icant gas emission detected at the inner disk by Spitzer 
(Najita et al. 2010). Thus, the depletion of dust optical 
depth appears to be orders of magnitude larger than any 
reduction in the gas accretion rate. This poses strong 
constraints on theories of the origin of pre/transitional 
disk structure, as we now show. 

3. METHODS: DISK-PLANET SIMULATIONS 

To simulate multiple giant planet systems, we use 
FARGO (Masset 2000). This is a two-dimensional hy- 
drodynamic code which utilizes a fixed grid in cylindrical 
polar coordinates {R, </>). FARGO uses finite differences 
to approximate derivatives, and the evolution equations 
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are divided into source and transport steps, similar to 
those of ZEUS (Stone et al. 1992). However, an orbital 
advection scheme has been incorporated which reduces 
the numerical diffusivity and significantly increases the 
allowable timestep as limited by the Courant-Friedrichs- 
Lewy (CFL) condition. Thus FARGO enables us to 
study the interaction between the disk and embedded 
planets over a full viscous timescale for disks whose inner 
radii are considerably smaller than the radial locations 
of the embedded planets, which is an essential require- 
ment for studying large gaps created by multiple planet 
systems. 

We assume a central star mass of IM0 and a fully 
viscous disk. We further assume a radial temperature 
distribution T = 221(i?/AC/)-i/2 K, which is roughly 
consistent with typical T Tauri disks in which irradiation 
from the central star dominates the disk temperature dis- 
tribution (e.g., D'Alcssio ct al. 2001). The disk is locally 
isothermal. The adopted radial temperature distribu- 
tion corresponds to an implicit ratio of disk scale height 
to cylindrical radius H/R = 0.029{R/AUf '^^. This dif- 
fers from the H/ R = constant assumption used in many 
previous simulations, which implies a temperature dis- 
tribution T oc R~^, which is inconsistent with observa- 
tions. Consequently, our assumed temperature distribu- 
tion with a constant viscosity parameter a {v = acl/il, 
where u is the kinematic viscosity, Cg is the sound speed, 
fl is the angular velocity) leads to a steady-disk surface 
density distribution E oc R^^ instead of the S oc R~^/'^ 
which would result from either assuming both H/R 

and M are constant, or both v and viscous torque (- 
2-KRYivR?dSl/ dR) are constant. This makes a significant 
difference in the innermost disk surface densities, and 
thus the implied inner disk optical depths in our models 
will be larger. 

We set a = 0.01 for the standard cases. Given our 
assumed disk temperature distribution, we set the initial 
disk surface density to be £==178 (R/AU)"^ (0.01/a) 
g cm^^ from i? ~ 1 — 200 AU, which yields a steady 
disk solution with an accretion rate M ~ 1O~^M0 yr~^, 
typical of T Tauri disks (GuUbring et al. 1998; Hartmann 
et al. 1998). The mass within the radius Rout is then 



M{R < Rout) = / 2TrRT,{R)dR 



With i?„i=l AU, i?ou/=200 AU, and a=0.01, the total 
disk mass is thus 0.025AfQ, reasonably consistent with 
observational estimates (e.g., Andrews & Williams 2005). 

We use logarithmic spacing for the radial grid to span 
the large dynamic range needed. The azimuthal reso- 
lution is 256 cells covering 27r radians, and the radial 
resolution is adjusted so that each grid cell is square 
{RAO^AR). 

Each planet in the disk can interact gravitationally 
with the other planets and with the disk. The inter- 
actions among multiple planets and the star are solved 
using a 5th order Runge-Kutta integrator. Each planet 



feels the acceleration from the disk given by 

S(r')(rp - rydS 

+ e2)3/2 ' 
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(2) 



whore the integral is performed over the whole disk area 
S. Here denotes the planet's position, and e is the 
smoothing length set to be 0.6H, where H is the disk 
scale height at the planet's orbit. 

To accurately simulate close encounters between plan- 
ets, an additional timestep constraint is applied. For 
each pair of planets wo calculate the minimum 

timestep 

1 /^27r , 

at = -r^mm { | , (3) 
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where = [G{mi + mj)/\ri — Tj]^]^^"^ is an estimate 
of the orbital frequency that each pair of planets (with 
masses m, and rrij) would have if they were in orbit 
around each other, with jr, — Vj \ being the absolute value 
of the distance between them. Throughout most of the 
simulations the time step size is determined by the nor- 
mal CFL condition, but when a close encounter occurs 
the interaction is resolved by allowing Eq. (3) to deter- 
mine the time step size. 

Accretion onto the planet is simulated by depleting the 
disk surface density within the planet's Hill radius, Rh = 
{Mp/3M^y/^rp, where Mp and arc the planet and 
stellar masses, respectively. This approach is required 
because of the necessity to evolve the disk-planet system 
for longer than the viscous time scale measured at the 
disk outer edge, which forces us to adopt a relatively 
low numerical resolution that prevents us from modeling 
accurately the gas flow onto the planet within the planet 
Hill sphere. Defining a dimensionless scaling parameter, 
/, we deplete the disk within 0A5Rh at a rate 
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and between 0.45 and 0.75 Rh at a reduced rate 
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Jffp 



(4) 



(5) 



where Tp is the orbital period of the planet^ . We term 
/ the planet accretion timescale parameter. During the 

intermediate stages of giant planet formation, after the 
formation of the solid core, but prior to the onset of rapid 
gas accretion, gas slowly settles onto the planet at a rate 
determined by the Kelvin-Helmholtz time scale. Depend- 
ing on the opacity of the planet envelope, this phase lasts 
for up to a few Myr (Pollack et al. 1996; Papaloizou & 
Nelson 2005; Movshovitz et al. 2010). Once the planet 
reaches a value of ^ 35 50 M0, rapid gas accretion 
onto the planet is able to ensue, relatively unimpeded 
by the thermodynamic evolution of the atmosphere, and 
the planet can grow quickly to become a giant, during 
this final phase of growth, however, the planet contracts 
rapidly, becoming much smaller than its Hill sphere, and 
further gas accretion must occur via flow through a cir- 
cumplanetary disk (Papaloizou & Nelson 2005). Simu- 
lations presented by Aylifi'e & Bate (2009) suggest that 

^ Note that our parametrization is different from the publicly 
available version of FARGO. 
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the circumplanctary disk has a radius Rep ~ Rh/^ and 
aspect ratio {H/R)cp — 0.5. Adopting these values, the 
viscous time scale for the circumplanetary disk 

Tacc ^ l/(57ra)rp = f/{2n), (6) 

from which we see that / = 2 /{5a). Varying the value of 
a between representative values of 10^^ and 10^'^, leads 
to variations of the accretion time of between ^ 6 and 
60 planet orbits. Since / is poorly constrained, we vary 
it by two orders of magnitude in this work: 1 < / < 100. 

Observations of exoplanets allow us to place some con- 
straints on the value of /. Once gap formation ensues, 
mass accretion into the planet Hill sphere is controlled 
by the rate at which material is supplied to the planet by 
viscous evolution of the protplanetary disc. A mass ac- 
cretion rate of 10~^ M0/yr implies growth of the planet 
to masses > 10 Mj over disc lifetimes of ~ few Myr if 
f — 1. Given that exoplanets with masses > 10 Mj are 
relatively rare, this suggests that / = 1 is a reasonable 
lower limit. A value of / = 100 implies mass accretion 
onto the planet at the rate ~ lO'^-lO'^o Mq/jt (lO'^- 
10~^ Mj/yr). Such low values of / can barely explain 
the numerous exoplanets with masses in the range 1-5 
Mj, suggesting that / = 100 is a reasonable upper limit. 

We initiate most of our simulations with relatively low 
planet masses (details are given later in this section), and 
the fact that planets grow slowly if / = 100 means that 
during our simulations rapid type I and type III migra- 
tion are likely to dominate the evolution of the planet 
orbits (Masset and Papaloizou 2003). To avoid this, we 
initiate simulations with / = 100 with the smaller value 
/ = 10 until a gap opens and the planet reaches 1 Mj, 
at which point / is switched back to 100. The mass 
and momentum of the depleted portion of the disk are 
added to the planet's mass and momentum, such that 
the planet's migration will be (modestly) affected by its 
accretion from the disk. 

As discussed by Crida, Morbidelli, & Masset (2007), 
a standard open inner boundary condition (Stone et al. 
1992) in a fixed 2D grid can produce an unphysically 
rapid depletion of material through the inner boundary 
in the presence of the planets. There are two reasons for 
this. First, due to waves excited by the planet, the gas 
in the disk can have periodic inward and outward radial 
velocities larger than the net viscous velocity of accret- 
ing material. Thus, with the normal open boundary, 
material can flow inward while there is no compensat- 
ing outflow allowed. Second, the orbit of the gas at the 
inner boundary is not circular due to the gravitational 
potential of the planets; again, as material cannot pass 
back out through the inner boundary, rapid depletion of 
the inner disk material is enhanced. As we arc inter- 
ested in the amount of gas depiction in the disk inward 
of the planet-induced gap(s) over substantial evolution- 
ary timescales, it is important to avoid or minimize this 
unphysical mass depletion. 

Crida et al. (2007) were able to ameliorate this prob- 
lem by surrounding the 2D grid by extended ID grids 
(see their Figure 5). We follow Pierens & Nelson (2008), 
who found reasonable agreement with the Crida et al. 
results while using a 2D grid only by limiting the inflow 
velocities at the inner boundary to be no more than a fac- 
tor (3 larger than the viscous radial velocity in a steady 
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where Vin and Rin arc the viscosity and radius at the 
inner boundary. Pierens & Nelson (2008) found satis- 
factory behavior for /5 = 5 in their simulations, based 
on a comparison with 1D-2D calculations. As our pa- 
rameters differ somewhat from those of Pierens & Nel- 
son, we also made a comparison with the 1D-2D setup 
in FARGO described by Crida et al. (2007), and found 
that (3 = 3 provided reasonable results. We adopted the 
normal open boundary condition as our outer boundary 
condition ^. The open outer boundary allows mass to 
leave the computational domain, and in strongly per- 
turbed disks containing massive planets this can cause 
a noticeable reduction in the disk mass over long time 
scales which in principle can reduce the gas accretion 
rate through the disk. But because this effect requires 
a long time to develop, its effect on our results is very 
modest. 

We do not use the 1D-2D grid method in this paper 
because in some cases our planet masses grow consider- 
ably larger than one Jupiter mass, and the non-circular 
motion of the gas close to the 2D boundary is significant 
enough so that an artificial gap is opened at the interface 
between the ID and 2D grids, which causes the code to 
crash. 

Since our numerical inner boundary is far away from 
the real disk inner boundary, we assign the ghost zone 
density to be 

S, = Si|i (8) 

Rg 

to simulate the power law section of the constant a disk 
similarity solution, where subscript g denotes the ghost 
zone and 1 the first active zone. 

In order to reduce problems related to relaxation from 
initial conditions, and to mimic the effect of accreting 
gas onto a protoplanet at the beginning of the rapid gas 
accretion phased discussed earlier in this section, we per- 
formed the following steps. First, we allowed for the 
planet mass to grow linearly from zero to its initial mass 
during the first 5000 years. For simulations where the 
planet accretes gas from the disc the initial mass is 0.1 
Mj, and is 1 Mj otherwise. During this period the planet 
does not feel the gravitational force from the disk and 
does not accrete gas. This step is necessary; otherwise 
the sudden presence of a one Jupiter mass planet in a 
Keplerian disk disturbs the disk so much that the ex- 
cited acoustic waves can be present in the disk for a long 
time ^, and these waves gradually deplete the disk mass 
by creating fiows through the outcir boundary. Then we 
allow the planet to start accreting at lO'^ yr. Finally, we 
allow the planet to feel the disk's gravitational force grad- 
ually, starting at 10^ yr and with full force at t = 2 x 10"* 
yr. This gradual ramping up of the gravitational effects 
is also important, since a sudden turn-on of gravity can 
kick the planet outside the gap it just opened. As well 
as being convenient from a computational point of view, 

^ This is different from the publicly available version of FARGO. 
^ Waves can reflect back from the boundary even with the open 
boundary condition. 
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this approach also ensures that gap formation arises both 
because of tidal torques and accretion onto the planet. 

Problems with boundary conditions also led us to per- 
form simulations without the indirect terms in the grav- 
itational potential experienced by the disk and planets. 
These terms arise when working in the non-inertial ref- 
erence frame; bascid on the central star, as they account 
for the acceleration of the central star by the disk and 
planets. The resulting non-circular motions can cause 
problems at both inner and outer boundaries of our cir- 
cular grid similar to those discussed above. Given the 
importance of these issues to our investigation of the disk 
surface density distribution exterior to gaps, we therefore 
did not include the indirect terms initially. We subse- 
quently recomputed the models P4AN and P4A10 (see 
following section) to include the indirect terms to en- 
sure that the planetary stability properties and the gap 
structures that we report in this paper remain effectively 
unchanged: wc do indeed find that this is the case. (The 
result for the P4A10 case with indirect terms is shown in 
Fig. 8) In performing these tests with the indirect terms 
we found it necessary to adopt the wave-damping bound- 
ary condition described in De Val-Borro et al. (2006) in 
order for the disk to remain well behaved at the out- 
side edge. Adopting this boundary condition has very 
little influence, however, on the structure of the gap or 
the stability properties of the planetary systems, so we 
believe that the results presented in this paper provide 
an accurate description of the features most important 
for comparing with observed pre/transitional disk struc- 
tures. 

4. RESULTS 

We performed simulations with 1, 2, 3 and 4 planets 
in the disk to study planet gap opening and gas flowing 
through gaps. The disk simulations evolve over 1 million 
years (Myr), which is the characteristic viscous timescale 
at i? ~ 100 AU, and is also shorter than the time scale 
over which most T Tauri disks are observed to clear. The 
properties of the simulations are summarized in Table 2. 

4.1. Viscous disks with one planet 

For comparison and test purposes, we first calculated 
a case without planets (case PN). The disk evolves as 
expected for a viscous disk with constant a, maintaining 
a surface density S oc and a nearly constant mass 
accretion rate of 1O~^M0 yr"^. 

We then placed a single IMj (ramping up to IMj dur- 
ing the first 5000 yrs as discussed in §3), non- accreting 
planet in the disk at 20 AU initially (case PIAN). Af- 
ter 0.5 Myr the planet has migrated in to ~ 8 AU (Fig. 
1). The azimuthally-averaged surface density within the 
gap is reduced by about an order of magnitude (solid 
curve in the left panel of Fig. 1). The resulting gap 
is relatively shallow, as expected due to the high disk 
viscosity (Lin & Papaloizou 1993; Crida & Morbidelli 
2007) and relatively large disk thickness at 20 AU where 
H/R ~ 0.061. The viscous criterion for gap formation 
is q > 40/ Re, where q is the planet-star mass ratio and 
Re is the Reynolds number. At 20 AU in a disk with 
a = 0.01 clear gap formation is only expected for a planet 
with mass Mp > l.bMj. The thermal condition for gap 
formation, which determines if the disk response to the 



planet is non linear, requires Rh > H at the planet lo- 
cation, which is only just satisfied for a IMj planet at 
20 AU. 

Because the planet is not allowed to accrete from the 
disk, and the tidal torques that it exerts are not strong 
enough to truncate the mass flow through the disk, ac- 
cretion flow is channeled through the gap and arrives 
in the inner disk. Due to the continuous replenishment 
from the outer disk beyond the planet, the inner disk 
can maintain the lO~*M0yr~^ accretion rate onto the 
central star (right panel of Fig. 1). As expected, the 
inner disk surface density is not strongly affected by the 
presence of the planet and is close to that obtained in the 
unperturbed viscous disk model PN described above. 

We next allowed the planet to accrete, at a rate con- 
trolled by the parameter / as described in §2. The planet 
starts out with a mass of O.lMj and is again placed at 
20 AU. If / = 1 (case PlAl), after the initial phase 
when the planet accretes material aroimcl it and opens 
a gap, the inner disk surface density decreases by one 
order of magnitude, and consequently the accretion rate 
onto the central star also decreases by one order of mag- 
nitude (lO^^^Mo yr^^) relative to the non-accreting case. 
This case is nearly equivalent to that studied by Lubow 
& D'Angelo (2006) ^, who similarly found that - 90% of 
the disk accretion flows onto the planet, leaving only ^ 
10% to pass through into the inner disk. On the other 
hand we have much higher disk surface densities inside 
the gap than do Varniere et al. (2006). As discussed in 
§2, this difference is almost certainly due to the differ- 
ence in the inner boundary condition. In test calculations 
performed which did not use the /3-limiter on the radial 
velocity (§2), we observed a similar rapid depletion of the 
inner disk (see also Crida & Morbidelli 2007). 

As the planet mass grows with time, the gap becomes 
wider and deeper due to the removal of gas from the disk 
and the increasing effectiveness of the tidal torques. As 
90% of the accreting disk mass is added to the planet at 
a rate ~ lO~^M0yr~^, the planet grows to lOMj after 
1 Myr (left panel in Fig. 2), so that the gap is much 
deeper than the case where the planet does not accrete 
(left panel in Fig. 1). 

If / = 10 (case PIAIO), the inner disk surface density 
and mass accretion rate onto the central star is reduced 
to 30% of the non-accreting planet case. Because the 
planet accretes more slowly than in the / = 1 case, the 
planet is less massive (right panel in Fig. 2) and so the 
gap is shallower. The eccentricity for this case is lower 
than for PlAl case (this is also true for multiple planet 
cases, as discussed below). 

If / = 100 (case PIAIOO), the inner disk surface den- 
sity and accretion rate onto the central star is barely 
affected by the planet (dotted curve in Fig. 1). Only 
10% of the mass flowing through the gap accretes onto 
the planet. Thus, it takes 1 Myr for the planet to grow 
to IMj. The accretion rates show an initial rapid decay 
due to the inner disk adjusting to the reduced flow from 
the outer disk; thereafter, the accretion rates show a slow 
decay or nearly-constant behavior, depending upon how 
fast the planet mass grows (which affects the mass trans- 
fer through the gap). With a as large as 0.01, the gap is 

® Although their / = 0.1, the radius within which mass is added 
to the planet is 0.2 Hill radii. 
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shallow and the planet is not fully locked to the viscous 
evolution of outer edge of the gap, as there remains sig- 
nificant mass flow through the gap. The resulting planet 
migration timescale is ~ 10^ yr (Fig. 2) while the viscous 
timescale at that radius is f .8 x 10'"' yr; thus the planet 
does not migrate very far during the viscous timescale, 
allowing the inner disk mass accretion rate to approach a 
steady state. The reason for the migration being so much 
slower than the viscous evolution rate is that the flow of 
gas through the gap causes a positive corotation torque 
to be exerted on the planet (Masset 2002), reducing the 
rate of inward migration. This timescale estimate also 
suggests that if a is as large as 0.01, the supply of gas to 
accrete onto the star in transitional and pre-transitional 
disks must come from the outer disk to be sustained for 
~ 1 Myr. 

4.2. Viscous disks with two planets 

Figure 3 shows the cases (P2A) with two accreting 
planets originally placed at 12.5 and 20 AU. The gap 
opened by the joint action of the two planets is consid- 
erably widcir than in the single planet case (by roughly a 
factor of three between the outer and inner edge of the 
gap). 

If neither of the planets are accreting (P2AN), the 
inner disk surface density (solid curve in Fig. 3) and 
the mass accretion rate onto the star are close to the 
no planet case as expected after a period of initial re- 
laxation. However, if both of the planets are accreting 
rapidly (/ = 1, P2A1) the accretion rate onto the cen- 
tral star and the surface density in the disk interior to the 
gap arc both reduced by almost two orders of magnitude 
(long-dashed curve in Fig. 3), as expected since each 
planet accretes roughly 90% of the viscous flow, leaving 
only 1% to make it through the gap. 

The planets' orbital properties for the P2A1 case are 
shown in the left panel of Figure 4. The two planets 
migrate inwards at a significantly slower rate than in the 
one planet case (Kley 2000; compare Fig. 2). Due to the 
common gap formed by two planets, the torque from the 
outer disk on the inner planet is almost zero. Although 
the outer planet tries to migrate toward the inner planet, 
the latter pushes the outer planet outwards by locking it 
into 2:1 resonance. The eccentricity of the planets is also 
driven up as they move into 2:1 resonance (Fig. 4), as 
seen in other simulations (e.g., Snellgrove et al. 2001; 
Kley et al. 2005; Pierens & Nelson 2008). The outer 
planet's mass grows faster than the inner planet because 
the outer planet is continuously being fed by material 
from the outer disk. 

With median and slow planet accretion rates (/=10 
and 100, cases P2A10, P2A100), the disk accretion rates 
onto the central star and the inner disk surface densities 
are reduced by factors similar to that of the equivalent 
single planet cases. Unlike the P2A1 case, both of these 
planets almost grow at the same rate (right panel of Fig. 
4) because there is enough mass flow passing the outer 
planet to feed the inner planet. 

4.3. Viscous disks with three and four planets 

Including more giant planets in the disk produces 
wider gaps, as expected (see Figs. 5 and 7 for the case 
with three and four planets, respectively). However, the 



greater the number of planets that are present, the more 
the accretion flow reaching the inner disk is reduced. Fast 
planetary accretion (/=1) quickly results in decreasing 
both the inner disk surface density and the mass accre- 
tion rate by more than two orders of magnitude, while 
median and slow planet accretion (/=10 and 100) lead 
to less than one order of magnitude reduction initially. 
There is a longer-term decay of the mass accretion rate, 
which is due to the fact that the inclusion of more plan- 
ets disturbs the disk more, and over long time scales this 
causes mass to be depleted from the disk as it escapes 
through the open outer boundary. 

The planets' orbital properties are shown in Figs. 6 
and 8. In most cases, the planets are trapped into 2:1 
resonances between each pair of adjacent planets. The 
middle panel of Fig. 8 displays the time evolution of 
the resonant angle (j)=2Xo-Xi-'cui associated with the 2:1 
resonance, where A/(Ao) and wi{xuo) are respectively 
the mean longitude and longitude of pericenter of the 
inner (outer) planet of the pair. Whenever the reso- 
nant angle librates between the interval [-7r,7r] the planets 
are trapped in a 2:1 resonance. The smaller the ampli- 
tude of libration, the deeper resonance locking. As in 
the two planet case (Figure 4), the planets' migration 
slows down due to resonance locking, possibly as a result 
of the growth of planetary eccentricity which can reduce 
the migration torques, and also due to the larger inertia 
associated with a resonant multi-planet system. 

There are two cases, however, where planets experi- 
ence a close encounter and are gravitationally scattered. 
The first case is P3A1 (left panel in Fig. 6), where the 
eccentricity of the middle planet is driven so high that 
its orbit overlaps with the inner planet's orbit, leading 
the inner planet to be scattered out from this system 
(in a hyperbolic orbit). The second case is P4AN (left 
panel in Fig. 8), where the innermost two planets slip 
out of 2:1 resonance leading to scattering. In this case 
the planets remain bound but switch their positions. It 
is clear that the disk model with a = 0.01 leads to cir- 
cumstances where three or four planet resonances can 
be maintained stably over Myr time scales, but can also 
lead to the break up of the resonances. This diversity 
of behavior appears to have a strong stochastic compo- 
nent, such that the stability of the planet system depends 
on the detailed history of disc-planet and planet-planet 
interactions. 

We also performed some limited calculations of three 

or four planets in disks with a lower viscosity (q;=0.002) 
and a higher disk mass to give the same disk accretion 
rate. In contrast to the lower mass a = O.Ol disks, these 
models all produced unstable planet systems which re- 
sulted in scattering and ejection of at least one planet 
from the system. This apparently occurs for two rea- 
sons. The first is that the low a disks result in deeper 
gaps, leading to a reduction in the eccentricity damping 
provided by the disk, thus favoring planetary instability. 
Secondly, the significantly larger disk mass probably has 
a more important dynamical influence on the resonant 
planetary systems, possibly helping to render the reso- 
nant configuration unstable. We will consider the effect 
of varying the viscosity in more detail in a future paper. 

5. DISCUSSION 
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Returning to the questions posed at the end of the 
Introduction, our findings are as follows: 

1) As seen in previous simulations, a single planet 
opens up a small gap. Multiple planets can open wider 
common gaps. However, in order to explain the prc- 
transitional disk gaps spanning almost an order of mag- 
nitude in radius without dust depletion (§2), we need as 
many as three to four giant planets, given the tendency 
of the viscously evolving disk to drive the planets into 
2:1 resonances. This is a large number of such planets, 
given current exoplanet statistics which do not include 
examples of systems containing more than two planets 
in resonance. One possibility is that resonant planet sys- 
tems which contain more than two planets are able to 
remain stable over significant time periods (~ 1 Myr) in 
the presence of the gas disk, whose contribution to eccen- 
tricity damping helps maintain dynamical stability. Once 
the gas disk dissipates, however, these resonant planetary 
systems may become unstable, leading eventually to ejec- 
tion of some members and the formation of more sparsely 
populated systems of planets on eccentric orbits, similar 
to those observed. According to this hypothesis, large 
systems of planets help explain the presence of the large 
disk gaps in pre-transitional disks, and also explain the 
population of eccentric giant cxoplanets. 

In principle, the gaps could be widened for fewer plan- 
ets if the orbits were eccentric; however, with large 
a=0.01, we do not see high eccentricities in stable sys- 
tems. 

2) A high accretion rate past the planetary gap can 
be maintained if accretion onto the planets is sufficiently 
slow. Planetary accretion rates of / = 10 still permit 
substantial accretion (>10% accretion rate outside the 
gap) past the planets to the inner disk. However, the 
reduction in surface density interior to the gap is in pro- 
portion to the reduction in mass accretion rate. This 
differs from the results of Varniere et al. (2006) but is 
consistent with the results of Crida et al. (2007), and is 
almost certainly the result of the improved treatment of 
the inner boundary condition. 

3) Carefully comparing our simulations with observa- 
tions indicates that the pre/transitional disk systems re- 
quire substantial reduction in dust opacities within the 
gaps (and inner disk if present) since the reductions in 
surface density there are not sufficient to explain the very 
small amount of dust required to fit the observations. 
For example, if we adopt the ISM opacity at 10 /xm (the 
assumption used by Espaillat et al. (2010) to estimate 
optically-thin dust masses) of k ~ 10 cm^ g~^, even the 
gaps with surface densities S > O.lgcm"^ are still op- 
tically thick, which is the case for all but the / = 1 
simulations. / = 1 simulations still have optical depth 
^0.1 and they are the cases with an extreme upper limit 
on the expected planetary accretion rate, which leads to 
low accretion rate onto the star (esp. with multiple plan- 
ets). Furthermore, the dust depletion must be larger in 
the inner disk than in outer regions. This follows from 
our adoption of a realistic temperature distribution and 
improved inner boundary condition, which results in a 
significantly higher surface density at small radii than 
at large radii; thus the inner disk will not be optically 
thin unless the outer disk is also optically thin, which 
disagrees strongly with transitional disk observations. 



One alternative is to place planets sufficiently close to 
the central star that the gap essentially extends either up 
to the dust destruction radius, or near to it; however, in 
that case it is difficult to extend the gap all the way out 
to 20 - 70 AU as observed (Table 1) when using only four 
planets. In principle it may be possible to accommodate 
additional planets at smaller orbital radii which truncate 
the inner disk to a smaller radius. The other alternative 
is to allow the planets to accrete more mass, depleting 
the inner disk; but then the mass accretion rate onto the 
central star is reduced by unacceptably large values. 

Protoplanetary disks are probably not fully viscous 
(e.g., Gammie 1996), but the existence of "dead zones" 
and/ or rediiced viscosity regions is unlikely to improve 
the situation, as they will lead to higher surface densities 
for a given mass accretion rate, and may reduce the sta- 
bility of systems of more than two planets in resonance. 

5.1. Gap structure and dust depletion 

A multi-giant planet model coupled with substantial 
dust depletion in the gap and inner disk does have an at- 
tractive feature. As pointed out in §2, there is some evi- 
dence - both from SED modeling and from ncar-infrarcd 
interferomctry - that the pre-transitional and transitional 
disks exhibit structure within the optically-thin disk gap 
or hole, such that the innermost gap region is much less 
cleared of small dust than the outer gap region. Schemat- 
ically, we might identify the highly-cleared region with 
the true planet-driven gap, and the inner, optically-thin 
dust region with accreting gas that has been strongly 
depleted in dust but not so depleted in gas surface den- 
sity. The idea is illustrated in Fig.9, which compares the 
inferred structure of one of the pre-transitional disk sys- 
tems with the azimuthally-averaged dust surface density 
for the P4A10 case. 

In Fig. 9, the disk optical depth at 10 /im can be 
estimated by multiplying the dust surface density (dust 
refers to small dust particles which contribute to 10 iim 
opacity) by 1000 cm^/g (which is estimated using the 
ISM opacity: a factor of 100 comes from the gas-to- 
dust mass ratio and 10 cm^/g is the ISM opacity at 10 
fixn). Thus, disk regions with dust surface density >10~^ 
g/cw? are optically thick. The dotted curve represents 
the dust surface density resulting from the model P4A10 
assuming there is no dust depletion, which is obtained 
by dividing the gas surface density from the P4A10 sim- 
ulation by the nominal gas-to-dust mass ratio ^100. In 
this case even the gap is optically thick. The solid curve 
represents the case where we notionally deplete the dust 
content inside the gap by a factor of 100. This can ex- 
plain the pre-transitional disks whose inner region within 
~ 1 AU changes from being optically thick to optically 
thin, with the region beyond being almost dust free. This 
comparison with pre-transitional disks set a lower limit 
on the dust depletion factor, while the transitional disks 
can set a higher limit. If we assume all the dust inside the 
transitional disk gap has been detected (which is reason- 
able since transitional disks are optically thin), we can 
estimate the dust surface density by solving 

I j.^^^^27ri?rfi? = Md , (9) 

Jo 

where Rout, .thin and M^; are given in Table 1, and is as- 
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sumcd to be Ed(AU)(R/AU)-i. The derived for GM 
Aur (<1 AU) using Table 1 is shown as the dashed curve 
in Fig. 9. By comparing with the dotted curve, we find 
the dust needs to deplete by a factor of ~10^. Thus we 
estimate the dust-to-gas mass ratio in the inner disk is 
between 10"^ and 10"^ of the ISM value (0.01) for our 
model to be consistent with pre/transitional disks. 

More broadly, a gap-opening perturbing body or 
bodies still remains a plausible explanation of the 
pre/transitional disks. This type of model naturally pre- 
dicts a sharp transition in surface density at the outer 
gap edge, which is needed to explain the strong mid- and 
far-infrared emission of these systems. While radially- 
dependent dust depletion might also mimic this effect, 
most T Tauri disks do not show this behavior. There are 
a significant number of T Tauri systems in which the disk 
is optically-thick at a wide range of radii, but which are 
much more geometrically flat based on their SEDs (e.g., 
the "group E" systems shown in Figure 7 of Furlan ct al. 
2006). These objects suggest that dust growth and set- 
tling can occur roughly simultaneously over a wide range 
of radii, without resulting in the abrupt change between 
optically-thin and optically-thick regions that occurs in 
the pre/transitional disks. In addition, the large outer 
disk gap radii (as much as 20-70 AU) pose a challenge 
for pure coagulation models at an age of ^ 1 — 2 Myr 
(Tanaka et al. 2005; DuUemond & Dominik 2005). 

The requirement of signiflcant dust depletion is not 
completely surprising, given the evidence for small dust 
depletion in many T Tauri disks without obvious gaps 
or holes (e.g., D'Alessio et al. 2001; Furlan et al. 2006). 
In addition, formation of giant planets via core accretion 
clearly requires dust growth. Indeed, an outstanding the- 
oretical problem has been to avoid clearing inner disks by 
ages of 1 Myr or less (e.g., DuUemond & Dominik 2005; 
Tanaka et al. 2005), resulting in suggestions that small 
dust must be replenished to some extent by fragmenta- 
tion as a result of collisions of larger particles. 

A hypothesis which explains many features of 
pre/transitional disks using both giant planets and dust 
depletion can be outlined as follows. Dust coagulation 
and growth in the disk interior to ^ 20 AU ensues, lead- 
ing to the formation of planetesimals, and eventually 
to the formation of a system of numerous giant plan- 
ets. This system of giant planets forms a large com- 
mon gap which covers radii from < 1 AU out to ~ 50 
AU. Remnant planetesimals within the gap region, whose 
collisions can act as a secondary source of dust, are dy- 
namically cleared out - preventing in situ secondary dust 
formation. After a viscous time scale corresponding to 
the size scale of the planetary system, the gas present in 
the planet-induced gap and the inner disk within 1 AU 
has originated largely from that part of the disk which 
lies out beyond the planetary system. This gas may be 
substantially depleted of small dust because of signif- 
icant grain growth at large radii, combined with filtra- 
tion of dust at the outer edge of the gap (Paardekooper & 
Mellema 2006; Rice et al. 2006). Furthermore, small dust 
particles may manage to pass through the gap, but they 
can quickly coagulate to big dust particles (DuUemond 
& Dominik 2005). In this scenario, gas which accretes 
through the system at late times is strongly depleted of 
dust, can sustain a significant accretion rate onto the cen- 
tral star, and can provide a dust-rich wall at the outer 



edge of the gap required by SED modeling. We will ex- 
plore this hypothesis in a forthcoming paper in which we 
include the effects of dust filtration. 

5.2. Num,erical limitations 

Our simulations are highly simplified, and this needs 
to be borne in mind when interpreting our results. The 
resolution we adopt is insufficient to resolve the gas flow 
within the planet Hill sphere, and so we are forced to 
adopt a simplifled approach to simulating gas accretion 
onto the planet, instead of simulating the accretion pro- 
cess directly. As such the detailed evolution of the gas 
flow in and around the planet Hill sphere may not be 
modeled with a high degree of fidelity in our simulations. 
The simulations are two dimensional, which may not be 
a good approximation early on when the planet masses 
are low, but improves as the planet masses become large 
and gap formation occurs. We have neglected a detailed 
treatment of the gas thermodynamics, adopting instead a 
locally isothermal equation of state. A proper treatment 
would allow the local disk temperature to be determined 
by a balance between viscous and stellar heating, and 
radiative cooling, and this might aff'ect details of the gap 
structure due to the changing optical depth of the gas 
and its thermal evohition there (although viscous heating 
is generally much smaller than stellar heating). Ayliffe 
& Bate (2010) suggest radiation due to circumplanetary 
disk accretion tends to suppress the spiral shocks and 
leads to a shallower gap (their Figure 18), which only 
increases the need for dust depletion within the gaps. 
Finally, we do not explicitly simulate the MHD turbu- 
lence that is believed to provide the effective viscous 
stresses in protoplanetary disks (Balbus & Hawley 1991). 
But it appears from previous simulations that adopting 
the usual a prescription gives results in broad agreement 
with MHD simulations when considering gap formation 
and giant planet migration (Nelson & Papaloizou 2003; 
Winters, Balbus & Hawley 2003; Papaloizou, Nelson & 
Snellgrove 2004). As such, we do not expect that the 
qualitative nature of our results have been compromised 
by the neglect of the above physical processes. 

6. SUMMARY 

In this paper we present 2-D hydrodynamic simula- 
tions to explore the possibility that the properties of 
pre/transitional disks are due to gap opening by plan- 
ets in viscous disks. With an improved inner boundary 
condition, we found that the surface density of the disk 
interior to the planet-formed gap is depleted by the same 
factor as the mass accretion rate through the gap and 
into the inner disk is reduced. Thus, the substantial ac- 
cretion rates of pre/transitional disks require substantial 
gas surface densities inside of planet-driven gaps. We 
also found that even multiple planets have difficulty in 
making large disk gaps and holes which are deep enough 
to explain the observations. Thus small grain depletion 
seems to be an essential part of the explanation of the 
structure of pre/transitional disks, and this may be ex- 
plained by the dust growth or the accretion of grain- 
depleted gas from large disk radii at late times and dust 
filtration. In addition, multiple planets might account for 
the inferred structure within the gaps of pre/transitional 
disks, in which the inner gaps are much less cleared of 
small dust than outer gap regions. 
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TABLE 1 
Stellar and Model Properties 
OF LkCa 15, UX Tau a, Rox 44, & GM AuR from Espaillat et 
AL. 2010 



LkCa 15 rX Tau A 
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(lAl Au 



M. (Mq) 
(K) 
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Optically Thick inner Disk 
<0.19 <0.21 <0.4 



^inner disk (AU) 
dinner disk (Mq) 



<2xl0- 



<6xlO- 



<8xl0" 



Optically Thick Outer Wall 

58 71 36 

Optically Thick Outer Disk 



R° „ (AU) 



20 



e 
a 



(Mq) 



0.001 
0.0005 
0.1 



0.001 
0.004 
0.04 



0.01 
0.006 
0.03 



0.5 
0.002 
0.16 





Optically Thin Dust 


in Gap/ Hole 




M (Mq) 


2x10"" 


2x10-" 


2xl0"i2 


<hnin (M^) 


0.005 


0.005 


0.005 


Q-max (MI^) 


0.25 


0.25 


0.25 
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4 


2 
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TABLE 2 
Models 



case name 




Rout 




Rp^ 


M p 


Acc 




AU 


AU 


g cm~'^ 


AU 


Mj 


f 


PN 


1.5 


200 


178 


- 


- 


- 


PlAN 


1.5 


200 


178 


20 


1 


No 


PIAIOO 


1.5 


200 


178 


20 


0.1 


100 


PIAIO 


1.5 


200 


178 


20 


0.1 


10 


PlAl 


1.5 


200 


178 


20 


0.1 


1 


P2AN 


0.75 


200 


178 


12.5/20 


1 


No 


P2A1UU 


0.75 


200 


178 


12.5/20 


0.1 


100 


P2A10 


0.75 


200 


178 


12.5/20 


0.1 


10 


P2A1 


0.75 


200 


178 


12.5/20 


0.1 


1 


P3AN 


0.5 


200 


178 


7.5/12.5/20 


1 


No 


P3A100 


0.5 


200 


178 


7.5/12.5/20 


0.1 


100 


P3A10 


0.5 


200 


178 


7.5/12.5/20 


0.1 


10 


P3A1 


0.5 


200 


178 


7.5/12.5/20 


0.1 


1 


P4AN 


0.5 


200 


178 


5/7.5/12.5/20 


1 


No 


P4A100 


0.5 


200 


178 


5/7.5/12.5/20 


0.1 


100 


P4A10 


0.5 


200 


178 


5/7.5/12.5/20 


0.1 


10 


P4A1 


0.5 


200 


178 


5/7.5/12.5/20 


0.1 


1 



^Thc initial planet radii. The radii can change with time due to migration and planet-planet interaction. 
^The initial planet mass. If Acc is not 'No', the mass can increase with time due to accretion from the disk. 
*^The accretion coefficient f is defined in Eqs 4 Sz 5. 
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Fig. 1. — Left: disk surfeice densities at 0.5 Myr for viscous disks (a=0.01) with one planet accreting at different rates: no planet accretion 
(solid curve), accretion on 100 (/=100) (dotted curve), 10 (/=10) (short-dashed curve), and 1 orbital timescales (/=1) (long-dashed curve; 
see text for further explanation). Right: the disk a<;cretion rates onto the central star for these four cases. A comparison between the left 
and right panels demonstrates that the depletion of the inner disk surface density results in a proportional decrease in accretion rate onto 
the central star. 
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Fig. 2. — 



The planet's semi-major axis, eccentricity , and mass as a function of time for the PlAl case (left) and PIAIO case (right). 



Transitional Disk Gaps and Planets 



13 




Fig. 3. — Left: disk surface densities at 0.5 Myr for two planets with no planet accretion (solid curve), /=100 (dotted curve), /=10 
(short-dashed curve), and /=! (long-dashed curve). Right: the disk accretion rates onto the central star for these four cases. 




Fig. 4. — The planets' semi-major axis, eccentricities, and masses as a function of time for the two axjcreting planets in the P2A1 case 
(left) and P2A10 case (right). The solid curves represent the outer planet while the dotted curves represent the inner planet. 
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Fig. 5. — Left: disk surface densities at 0.4 Myrs for three planets with no planet accretion (solid curve), /=100 (dotted curve), / = 10 
(short-dashed curve), and / = 1 (long-dashed curve). Right: the corresponding disk accretion rates. 




Fig. 6. — The planets' semi-major axis, eccentricities, and masses as a function of time for three accreting planets in P3A1 case (left) 
and P3A10 case (right). 
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Fig. 7. — Left: disk surface densities at 0.4 Myrs for four planets with no planet accretion (solid curve), /=100 (dotted curve), / = 10 
(short-dashed curve), and / = 1 (long-dashed curve). Right: the corresponding disk accretion rates. 
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Fig. 8. — The planets' semi-major axis, eccentricities, resonant angles (2:1), masses with time for four planets in P4AN case (left) and 
P4A10 case (right). In the panel of the right figure, black, red, and green curves represent the resonant angle between the pairs of 4/3, 
3/2, 2/1 planets (4 is the outermost one.). 
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Fig. 9. — The pre-transitional disk structure taken from Espaillat et al. (2010) compared with the azimuthal-averaged dust surface density 
of the P4A10 case. The dotted curve represents the case with no dust depletion and the solid curve represents the case with the dust at 
the inner disk depleted by a factor of 100. The solid curve can explain the pretransitional disk structure as indicated by the arrows. The 
dashed curve is the dust surface density from the transitional disk GM Aur estimated by Eq. 9, which requires the dust to deplete by a 
factor of 10^ by comparison with the dotted curve. 



